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Some Inversion Formulas for the Cone 

Transform 

Fatma Terzioglu* 


Abstract 

Several novel imaging applications have lead recently to a variety 
of Radon type transforms, where integration is done over a family 
of conical surfaces. We call them cone transforms (in 2D they are 
also called V-line or broken ray transforms). Most prominently, they 
are present in the so called Compton camera imaging that arises in 
medical diagnostics, astronomy, and lately in homeland security ap¬ 
plications. Several specific incarnations of the cone transform have 
been considered separately. In this paper, we address the most gen¬ 
eral (and overdetermined) cone transform, obtain integral relations 
between cone and Radon transforms in M”, and a variety of inversion 
formulas. In many applications (e.g., in homeland security), the signal 
to noise ratio is very low. So, if overdetermined data is collected (as 
in the case of Compton imaging), attempts to reduce the dimension¬ 
ality might lead to essential elimination of the signal. Thus, our main 
concentration is on obtaining formulas involving overdetermined data. 


1 Introduction 

In this paper, we study the so called cone transform, where a function on M”' 
is integrated over various conical surfaces (in 2D, the names V-line transform 
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and broken ray transform are also used). Such transforms arise in a variety 
of new imaging techniques, e.g. in optical imaging [B], but most prominently 
in the so called Compton camera imaging, which we will briefly explain now. 
The conventional gamma cameras used in medical SPECT(Single Photon 
Emission Tomography) imaging determine the direction of an incoming 7- 
photon by ’’collimating” the detector (see Fig. [Deleft)). This considerably 
decreases the efficiency, because only a small portion of the incoming 7- 
rays passes through the collimator [ 3 ]. Thus, the acquired signal is weak 
and statistically noisy. The situation is similar in astronomy and even more 
severe in homeland security applications [11121 [mi 26 ] . 

On the other hand, Compton cameras utilize Compton scattering (see Fig. 
[Upright)) and use electronic rather than mechanical collimation to provide 
simultaneous multiple views of the object and dramatic increase in sensitivity 

ra- 



Figure 1: Left: Collimation. Right: Compton Scattering. 


A Compton camera consists of two parallel detectors (see Fig. [ 2 ]). When 
the photon hits the hrst detector, where its position u and energy Ei are 
recorded, it undergoes Compton scattering. Then, it is absorbed in the 
second detector where its position v and energy E2 are again measured. The 
scattering angle and a unit vector ft are calculated from the data as follows 
(see e.g. [S]): 


cos if 


mc^Ei 

{El + E2)E2 


/S 


u — V 
u — V 


( 1 ) 


Here, m is the mass of the electron and c is the speed of light. 

From the knowledge of the scattering angle fj and the vector ft, we con¬ 
clude that the photon originated from the surface of the cone with central 
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axis / 9 , vertex u and opening angle 'ijj (see Fig. H]). Therefore, although the 
exact incoming direction of the detected particle is not available, one knows 
a surface cone of such possible directions. One can argue that the data pro¬ 
vided by Compton camera are integrals of the distribution of the radiation 
sources over conical surfaces having vertex at the detector. The operator 
that maps source intensity distribution function f{x) to its integrals over 
these cones is called the cone or Compton transform. The goal of Compton 
camera imaging is to recover source distribution from this data [1]. 



Figure 2: Schematic representation of a Compton camera. 


In the Compton camera imaging applications mentioned above, the vertex 
of the cone is located on the detector plane, while in other applications 
vertices are not restricted, although some other conditions are imposed on 
the cones. We thus hnd it useful to understand analytic properties of a 
more general cone transform, where no restriction on the vertex location is 
imposed. This is the transform addressed in this text with the hope that it 
can be useful for more restricted versions. As for instance Remark [ 15 ] shows, 
one indeed arrives at applications to the Compton imaging. 

The problem of inverting the cone transform is over-determined. For 
instance, the space of 2 D cones with vertices on a linear detector array is 
three-dimensional, and the space of 3 D cones with vertices on a detector 

^It is planned to address these applications in detail elsewhere. 
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surface is five-dimensional. Without the restriction on the vertex, the di¬ 
mensions are correspondingly four and six. One thus is tempted to restrict 
the set of cones, in order to get a non-over-determined problem. There exist 
several inversion formulas of this type (e.g. [31I1KZIEI]). However, as we have 
already mentioned, when the signals are weak (e.g. in homeland security ap¬ 
plications (e.g., PP), restricting the data would lead to essential elimination 
of the signal. We thus intend to use the full data set. 

Probably, the hrst known analytical reconstruction formula in 3D was 
given in [ 1 ], where the authors considered cones with vertical axis only. The 
papers [3l[T3] contain spherical harmonics expansion solutions. Another in¬ 
version formula for cone transforms on cones having fixed central axis and 
variable opening angle is provided in . The paper [ 21 ] presents two recon¬ 
struction methods for two Compton data models. The complete set of data 
was used in [1^116] . Inversion formulas for n-dimensional cone transform 
over vertical cones are provided in [9l|T0]. All these works only addressed the 
cones with the vertex on the detector. Inversion algorithms for various 2D 
cone transforms are given in [51[^I5] [T2] [TB]. 

In this paper, we derive various inversion formula^ for the full data cone 
transform in M”. In Section 2, we define the cone transform and state its 
basic properties. In Section 3, we obtain an integral relation between the cone 
and Radon transforms in M"' and deduce from it an inversion formula for the 
cone transform. In Section 4, we provide a different inversion formula derived 
from another integral relation between the cone and Radon transforms in 
M"". Both of these formulas provide reconstructions only at vertices of the 
cones which is an inconvenience for Compton imaging. However, the integral 
relation provided in Section 4 also enables us to associate the cone transform 
with the cosine transform. This result is given in Section 5, and through 
this relation, we obtain the Radon transform explicitly in terms of the cone 
transform in Theorem 14 which leads to a variety of inversion algorithms 
from Compton data as discussed in Remark 15. The results of a numerical 
simulation for n = 2 are also provided. In Section 6, we investigate the 
relationship between the cone transform and spherical harmonics. Finally, 
we prove some auxiliary technical results in Section 7. 

^The reader should recall that it is common to have a variety of different inversion 
formulas for Radon type transforms, which are all the same for perfect data, but react 
differently to unavoidable errors in data 1141120] . Having such a variety is even more 
important when dealing with overdetermined data, as in Compton imaging. 
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2 Definition and Basic Properties of the Cone 
Transform 

A round cone in MT" can be parametrized by a tuple where u G M"' 

is the cone vertex, vector (3 G is directed along the cone’s central axis, 
and '0 £ (0)^) is the opening angle of the cone (see Fig. [2]). Then, a point 
X G lies on the cone iff 

{x — u) ■ (3 = \x — u\ cos'ip. (2) 

The n-dimensional cone transform C maps a function / into the set of 
its integrals over the circular cones in R”. Explicitly, 

Cf{u,(3,i))= j f{x)dx (3) 

{x—u)-l3=\x—u\ cos -i/; 

where dx is the surface measure on the cone. 

The n-dimensional vertical cone transform maps a function / into the 
set of its integrals over the cones having central axis parallel to the a:,i-axis, 
and thus the vector (3 is equal to e„ = (0,..., 0,1) G R”. It can be written in 
terms of the spherical coordinates. Namely, 

OO 

C/(M,e„,'0) = j j f{u + p{{sm.ip)u:,cosip)){psm.f)Y~‘^dudp. (4) 
0 

In two dimensions, the equation ([2]) describes two rays with a common 
vertex (see Fig. [3]). A cone in two dimensions can be parametrized by a point 
M G R^ that serves as its vertex, an opening angle tjj G (0, tt) and a vector 
(3 = (3{(j)) = (sin 0, cos 0) G directed along the central axis. 
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Figure 3: A Cone in 2-diniensions. 


Then, the 2D cone transform of a fnnction / G is given by 

OO 

Cf{u, 13, 'll:) = C/(m, /?(0), -0) = j f{u + r(sin(r/> + 0), cos(0 + (j))))dr 

0 

OO 

+ J /(m + r(—sin('0 — 0), cos('0 — 0)))(ir. 

0 

( 5 ) 

As a straightforward calculation shows, analogously to the Radon trans¬ 
form, cone transform has an evenness property, and is shift and rotation 
invariant: 

Lemma 1. Let f G u G M"', 0 G 5"“^ and ifj G (0,7r). Then, 

(^) 

C/(m,- 0,0) = C'/(M,0,7r-0). (6) 

(ii) Let Ta be the translation operator in R"", defined as Taf{x) = f{x -b a) 
for a G R". We define 

:= Cf{u + a,/3,fi). 

Then, 

TaC = CTa. 
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(in) Let A be an n x n rotation matrix and M^fi^x) = f{Ax) be the corre¬ 
sponding rotation operator. We define 


:= Cf{Au, A/3,'ll:). 


Then 


MaC = CM A. 


3 Inversion of the Cone Transform 

In the following, we investigate the relation between the cone and Radon 
transforms and provide various analytical inversion formulas for the n-dimensional 
cone transform. 

We hrst recall that the n-dimensional Radon transform R maps a function 
/ on into the set of its integrals over the hyperplanes of M”. Namely, if 
u G S"'~^ and s G M, 



( 7 ) 


In this setting, the Radon transform of / is the integral of / over the hyper¬ 
plane orthogonal to u with signed distance s from the origin. 

The Radon transform is invertible on namely 


/ = a<n. 


( 8 ) 


Here, R* is the back projection operator, and a < n, is the Riesz potential 
acting on a function f{u) as 


(7“/)(0 = ie|-“/(0, 


where / is the Fourier transform of /. For instance, when n is odd, C ” is 
simply the differential operator 


jl-n _ 


with A being the Laplacian (see e.g. [2U]b 
Theorem 2. Let f G iS(R"'). Then, 
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(i) For any u eMT' and (3 E ^, we have 

r(^) 


Cf{u,l3,'ip)dijj = 


r(^) 


2yj-(?l i)/2 


Rf{u,u-u)du = /?#/?/(«). 


gn-l 


( 9 ) 


(ii) Let a function /i : S'” —)■ M &e such that J fi{(3)d(3 = 1. For any 

gn-l 


f E S(R" 


/(«) = 


7r-”/2r(f; 


2r(n- 1) 



ji '^(jf{u,/3,ip)ia{/3)d'iljd/3. (10) 


S"-i 0 


Remark 3. 

(i) One notices that according to the inversion formula IflU) consists 
of a backprojecting of the cone data, followed by a filtration (i.e., is 
what is called a FBP type formula). 


(ii) One can choose fi{/3) to be egual to a delta-function, which would elim¬ 
inate integration with respect to f in m- However, if the signal is 
very week, eliminating almost all values of (3 would lead to elimination 
of the signal. Thus weighted integration with respect to (3 allows for 
accounting for all data collected. 


Proof. We first prove the theorem for dimensions n > 3. 

TT TT OO 


j Cf{u,en,f)df) 

0 



f{u + p((sin'0)a;, cos'0))(psin'0)” ^doodpdf 


0 S’! 


= / / f{u + p(T)p'^ ^dpda = j f{u + x)\x\ ^dx =-^^^^R*Rf{u), 

5^-1 0 R" 

The last eqnality is dne to m Chapter 2, Theorem 1.5] (see also Corollary 
[22]). As both R and R"^ commute with rigid motions in R”, we obtain for 
any (3 E S”“^, 

TT 

j Cf{u,l3,'ip)di) = j^^^R*Rf{u). 


0 







Thus, for any function /x on S'"" ^ such that / fi{f3)df3 = 1, we have 


Sn-l 


r r 1 rf—i 

J j Cf{u,(l,i,)md'<pdp = -—R*Rf(u) = 

S"-i 0 


Note that the last equality follows from the area formula for the n-sphere, 
that is 


1 ^ 


n—1 


27r"/2 

w 


( 11 ) 


Using ([8]) with a = n — 1, and utilizing the duplication formula (see e.g. ra) 


r(5)r(2 +1) = 2'-='v^r(2z), 


( 12 ) 


we conclude that 


/(«) = 


2 r(n - 1 ) 



r-^Cf{u,/3,ij)ij{/3)dijdl3. 


0 


For the 2-dimensional case, we only need to provide the proof of ijOl) , 
since the rest of the proof stays the same. Assume for now that u = 0. By 
definition of the 2D cone transform, we have 


j Cf{D,l3{(f)),'ip)di) 
0 



0 0 
TT OO 



/(r sin('^ -1- 0), r cos('0 -|- (j)))drdip 
/(—r sin('0 — 0), r cos(0 — 0))dr(i0. 


0 0 


Changing variables, we obtain 


TT 

J f{r sin(0-f 0), r cos(0-1-0))(i0 
0 


71+0 

//(rsm^.rcos^)#, 

0 
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and 


Thus, 


/(—r sin('0 — 0), r cos('0 — = J /(rsin'0,r cos'0)(i'0. 

—7r+0 

TT OO 7r+0 

J Cf{0, (3{(j)),^p)d^p = J J f [r sm. cos'ip)d'ipdr. 

0 0 —7r+0 

TT 

Changing variables by letting 6 = — — 'ip and using 27r-periodicity of sine and 
cosine functions, we get 

n+rf> 2 tt 


J f {r sin'IpiT COS Ip)dij) = J f(r cos 9, r sin 9)d9 = J f{r cos 9, r sin 9)d9. 

-Tr+<p 0 

Therefore, 


27r OO 



1 


2it 


CI3{(p),'ip)d'ip = I I f{r cos 9, r sin 9)drd9 = - / Rf{9,0)d9, 


0 0 


where the last equality follows by letting n = 2 and p = 0 in fl2^ . Now, 
using the shift invariance of both cone and Radon transforms, we conclude 
that 

TT 

j Cf{u,/3,'ip)d'ip = ^ j Rf{u,u-u)du = ^R*Rf{u), 

0 51 

which is ([9]) with n = 2, so we are done. □ 

Corollary 4. For n = 3, the formula (na becomes 



0 

where A acts on the variable u. 


AC f{u, 13, ip)ii{l3)d'ipdl3, 
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4 An Alternative Inversion Formula 


For the derivation of an alternative inversion formula, we need the following 
relation between the cone and Radon transforms. 

Theorem 5. Let f G For any n G M"' and (3 G , we have 

TT 

j (g/(M,/j,^)sin^d^ = j Rf{oj,oj ■u)\oj ■ l3\doj, (13) 

0 

where |S'"'“^| denotes the area of the sphere . 

As in the case of the Radon transform, invariance properties play a key 
role in the inversion of the cone transform. In fact, due to rotational invari¬ 
ance, it suffices to prove flT^ only for the vertical cone transform. Moreover, 
shift invariance enables us to consider vertical cones having vertex at the 
origin only, that is m = 0 . 

Proposition 6. For any f G we have 


TT 

J (F/(0,en,^)sin^d^ = J Rf{u,0)\u-en\duj. (14) 

0 S "-1 

For the proof, see Section 17.21 


Proof of Theorem\^ We will use Proposition [6] and the properties of the 
cone transform to deduce Theorem [5l We hrst remind that the Radon trans¬ 
form commutes with shifts and rotations, that is R(Tuf){oj, s) = Rf{u;,s + 
oj ■ u) and MARf{oj, s) = Rf{Au, s) = R{MAf){oj, s). 

As cone transform also commutes with shifts. Proposition |6] implies that 


Cf{u,en,fj)simjjdxp= / C{Tuf){0, en,fj) simpd'ijj 



Sn-l 


R(Tuf){uj,0)\u ■ en\duj 



Sn-l 


Rf{u),u ■ u)\u ■ en\duj. 
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Next, for (3 E let A be the rotation matrix such that {3 = Aen and 

X = A~^u. As cone transform commutes with rotations, we further have 


J C f{u, I3,'ip)sinil)d'ip 
0 


j C{MAf){x,en,'ijj) simjjdxp 
0 



Sn-l 


R{MAf){oJ,oj ■ x)\oj ■ en\du. 


Due to the rotational invariance of the Radon transform, we have 


J R{MAf){oo,uj ■ x)\uj ■ en\duj = J MARf{oo,uj ■ x)\uj ■ en\duj 

Sn—1 1 

= J Rf{Au,uj ■ x)\uj ■ en\du = J Rf{Au,uj-A~^u)\uj-A~^(3\du 
— 1 — 1 

= J Rf{Auj,Auj-u)\Au-f3\du= J Rf{uj,u-u)\uj-f3\du, 

— 1 gn—1 


The last equality is due to the rotational invariance of the Lebesgue measure 
on the sphere. Hence, we obtain flT^ . □ 

Remark 7. As it will be mentioned in SectionlBi the assumption f G 
can be significantly weakened. The same applies to Theorem\B 

The equality flT^ enables us to invert the cone transform by utilizing the 
inversion formulas for the Radon transform. 

Theorem 8. Let f G For any u G M"', we have 


/(«) = 


r^(^) 

27r"-r(n) 



ji f{u, sin'fidfidfi. ( 15 ) 


S"-i 0 


Proof. Integrating both sides of (lT3ll with respect to (3 over S'” , we obtain 



C f{u, sin'fid'ipd[3 = 


TT 


5^1-1 0 




Rf{u,u-u) / \uj ■ (3\d(3du. 

5n-l 
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Using the rotation invariance of the Lebesgue measure on the sphere, for any 
uj G we compute 


\u: ■ f3\df3 = 


gn-l 

Thus, we get 



cos(/)|(sin0)"' ^d(j)d6 = 




n—2 1 


n — 1 


0 



Cf{u, (3, 'ip) sin'ipd'ipd(3 = 


71 


n—2 1 


5^-1 0 


n-l 


Rf{oj, oj ■ u)du 


f ^*mu) ^ -3P!pUH*mu). (16) 


n - 1 


5^-1 


2n-ir2Wy 


Note that, for the evaluation of the constant, we have used the area formula 
for the n-sphere, flTT]) and the duplication formula flT^ . Now, using formula 
dHj) with a = n — 1, we obtain the result. □ 

Corollary 9. For n = 3, the formula flT^ reads as 


fiu) = 


-1 

47|-3 



ACf{u, (3, tp) sin tpdpjd(3, 


0 

where A acts on the variable u. 


5 Relation of the Cone Transform with Co¬ 
sine Transform. Other Inversion Formulas 

The main goal of this section is to derive a formula which is applicable in 
Compton imaging. This is achieved in Theorem 14 and Remark 15. We 
start, however, with a relation between the cone transform and the cosine 
transform which is dehned as follows: 

Definition 10. The cosine transform of a function f G C{S'^ is defined 
by 

j f{(^)W (17) 
s~-i 

for alloje S^-\ 
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Now the relation fll3p can be written as 


€{R{Tuf)){P) = j R{TJ){uj,Q)\uj ■ (3\du j Cf{u,(5,'ilj)sui'iljd'ilj. 

5"-i 0 

( 18 ) 

The cosine transform is a continuous bijection of to itself 

(see e.g. [7], [22])- Since, for any / G Rf{u,0) is an even function 

in (7°°(5'’^“^), we can recover the function R{Tuf) by inverting the cosine 
transform. Before stating this inversion formula, we recall the dehnitions of 
the Beltrami-Laplace operator and the Funk transform. 

Definition 11. Let f G The Beltrami-Laplace operator As on 

is defined by 

(A5/)(n) = l^r(A/)(:^), (19) 

fI 

~ X 

where fix) = fi- —r) is the homogeneous extension of f to M", and A is the 
\x\ 

Laplace operator on M"". 

Definition 12. Funk transform of a function f G is defined by 

Ff{e) = J /(cr)de(T = J f{a)de(T. (20) 

{(T£S^-^-.d{(T,e)=-K/2} 

Here, d{a, 6) = arccos{a ■ 9) is the geodesic distance between the points a and 
9 in S^~^, and d^a stands for the 0{n)-invariant probability measure on the 
{n — 2)-dimensional sphere S"‘~^ fl 9^. 

Theorem 13. Let g = €f, f G C'“g„(S'”“^). Then, if n is odd, 

{ _27r(2-’^)/2 r 1 1 f 

Sn-l ) gn-l 

( 21 ) 

with r = (n + l)/2, and if n is even, 

7r2"“^ 

/ = cPr{As)Fg, c = --—, (22) 

1 (n — Ij 
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with r = n/2, where F is the Funk transform and 


r—1 

P,(As) = 4-^ n +{2k-l){n-l- 2k )], 

k=0 


with As being the Beltrami-Laplace operator on 5'"“^. 

Thus, we can find RTuf explicitly for all u G M”: 
Theorem 14. Let f G For any u G M" and u G 5"“^, 

(i) if n is odd, 

Rfioj, u ■ u) 


-27r-”/2 


P{n+l)/2{As) 



C f{u, (3, fj) log -j- — sin fjdipdfd 


1-1 0 


uj- /3\ 


r(!i±i: 

„(n+l)/2 



Cf{u,f3,ip) sin tjjdiljdft, 


(23) 


5^-1 0 


(a) if n is even, 


_on—1 r 

Rf{0J,0J- u) = _ 1) y Pn/ 2 {As)F{Cf){u,u,'il))sin'ipd'il), (24) 

0 

where F and Pr{As) are given as in Theorem [XM. and both of them act on 
the variable oj. 

Proof. The result follows by applying inverse cosine transform d^T]) and fl22p 
to equality flTS]) . □ 

Remark 15. 

(i) For any oj G 5"“^ and s G M, the Radon transform Rf{u,s) of a 
function f G S(ML) can be computed using formulas fl23l) and fl2Tl) . 
if for any {u,s) one has access to a cone vertex (= detector location) 
M G M”' such that u ■ u = s. For instance a line (curve) array of 
detectors should be sufficient. Thus, Theorem [j^l together with formula 
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(]5j) should provide inversion formulas for the cone transform that are 
applicable to Compton camera data. This idea leads to a variety of new 
inversion formulas for Compton camera imaging, which will be derived 
and applied elsewhere. 

(a) We applied this approach to some 2D examples. Figures\^and\^show 
the reconstructions of some phantoms from their projections collected by 
four Compton cameras placed along the sides of a sguare. We simulate 
analytically the Compton projection data of the phantoms and then use 
formula to convert them to Radon projections. Finally, the filtered 
back-projection is applied to invert the Radon transform and obtain the 
reconstructions. 



Figure 4: Left: The phantom is the characteristic function of a circle having density 
1 unit, radius 0.5 unit and centered at (0,0). Right: 256x256 image reconstructed 
from the simulated Compton data using 257 detectors per side and 200 counts for 
the angles ft and xf each (see Fig. ED. 
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Figure 5: Left: The phantom is the sum of the characteristic functions of two 
intersecting circles having densities 0.3 and 0.7 units, radii 0.5 and 0.3 units, 
and centered at (0,0) and (0.5,0). Right: 256x256 image reconstructed from 
the simulated Compton data using 257 detectors per side and 200 counts for the 
angles /? and '0 each (see Fig. El). 


6 Relation of the Cone Transform with Spher¬ 
ical Harmonics 


Utilizing the relation of the cosine transform with spherical harmonics, we 
can relate the coefficients of the spherical harmonics expansion of the cone 
and Radon transforms. 

Lemma 16. Let g G Then, 



Cf{u, = TT / Rf{u,uj-u)€g{uj)du. (25) 


0 


Proof. Multiplying both sides of flT^ with g{/3) and integrating with respect 
to /3 over S^~^, we have 



Cf {u, Id,'ip)g{l3) sin'ipdfjdld = 


TT 


0 




Rf{uj,uj-u) / g{f3)\uj ■ f3\d(3du 


= TT / Rf{u,uj-u)€g{u)du. 

Sr^-l 
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□ 


The spherical harmonics are known to be the eigenfunctions of the cosine 
transform. This follows from the Funk-Hecke Formula: 

Theorem 17 ( [H]). (Funk-Hecke Formula) Suppose f{t) is continuous for 
t G [—1,1]. Then, for every spherical harmonic of degree m andoj G 

j f{uj-a)Ym{a)da = X^Ym{uj), (26) 


with 

1 

A^ = |5"-'| j f{t)PUt){l-t^Y'^-^y^dt, 

-1 

where Pmit) are the Legendre polynomials, see fW) . 

Corollary 18. For every spherical harmonic Ym of degree m, m = 0, 1, 2, 
and every oj G 5"“^, 

€F^(a;) = (27) 

where Xm is given as in Funk-Hecke Formula for f(t) = |f|. 

Now, we can establish the following relation. 

Proposition 19. For every spherical harmonic Ym of degree m. 


TT 

J J Cf{u,/3,'ilj)Ym{/3) sin fjdfjd/3 
0 


TlXm 


Rf{u,u-u)Ym{uj)du. (28) 


5n-l 


In particular, for m = 0, we obtain (USD. 

Proof. Letting g = Ym in (l25D . and using (1271) . we get (1281) . Then, the 
equation (IT6i) follows from direct calculation. □ 

Remark 20. As the relation fl28ll gives the spherical harmonics coefficients 
of the function Rf{u, u-oj), one can recover it for all u ^ ML and oj G 
Then, any inversion formula for the Radon transform (jHD would reconstruct 
the function f. This can be considered as an analog of Cormack’s method 

f2Uf. 
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7 Proofs of Some Auxiliary Statements 

7.1 An Integral Relation for the Radon Transform 

Lemma 21. For any f G u G M”, and p G M, 


Rf{u,p + u-u)du = 15 


‘71 — 2 I 



f{u + ru){r‘^—p‘^Y"^ ^^^“^rdrduj. (29) 


gn-l 


S" 1 IpI 


Proof. Due to the shift invariance of the Radon transform, it suffices to prove 
the lemma for u = 0 only. Let F be the spherical mean-value of /, i.e., 


F{r) = 




f{ru)du. 


The rotational invariance of the Radon transform implies that it com¬ 
mutes with the spherical mean-value operator. Thus, 


F{p) := RF{u,p) = 




Rf{^,p)dC 


gn-l 


On the other hand, if {u,Ui ,..., is an orthonormal system in M”, 


OO OO 


F{p) = 


— OO —OO 

OO OO 


— OO —OO 


F[piJ + tiUJi + • • * + 

R (y\/p‘^+tl - dti...dtn-l 

as F is radial. Letting x = + • • ■ tn-iCOn-v have 

F dti...dtn-i = j Fi^^/p"^ + \x\‘^)dx 


OO OO 


— OO —OO 


= | 5 ”- 2 | / F{^p^ + R)r-^dt. 
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Finally, letting r = obtain 


y F{\Jp^ + = J F(r)(r^ — ^'^^'^rdr 

0 IpI 

CXD 

1 





/(ra;)(r^ — p^)*-” ^'^^‘^rdrduj. 


Hence, the result follows. 

Corollary 22 ( [20]). Letting p = 0 in (12^ . we obtain 


□ 


R*Rf{u)= J Rf{u,u-u)du=\S'^ J J ‘^drdu 

gn-l g-n-l 0 

= |5’^-2| [ f{u + x)\x\-^dx = \S^-^\{\x\-U f){u). (30) 


7.2 Proof of Proposition [6] 

We first prove the proposition for n = 2. By dehnition of the 2-dimensional 
cone transform (l5l), we have 

TT TT CO 

j C/(O,e 2 ,' 0 ) sm.'ipd'ip = j J f {r simjj, r cos ^jJ) simp dr d^p 

0 0 0 

TT OO 

+ yy /(—r sin'^, r cos V’) sm.'ipdrd'ip. 

0 0 

Changing variables by letting r ^ —r and "0 —)■ vr — -0, respectively, we 
obtain 

TT CO TT 0 

/(..inV-,rco.V-).ini.<i--# = j j /(-.=inV-, -r co.^) 

0 —CO 
TT 0 

f{—r sin 0, r cos (/>) sin (pdrdcp. 

0 —CO 



0 0 
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Therefore, 



C f{0, e 2 ,'i/j) simjjdijj = / / f{—r sin ip, r cosip) sin ip dr dip 


0 —oo 


= / Rf {uj {ip), Q) sin Ip dip 


where u:{ip) := (cos "0, sin'0). Now, the evenness property of the Radon 
transform implies that 


TT 


0 


Rf{uj{ip),R) sinipdip 


J Rf {u{ip + Tr),0) sinipdip 
0 


Hence, we get 


2tt 

J Rf{u{(p),0) sin(pd(p. 


TT 


TT 27r 

j Cf{0,e2,ip) sinipdip = ^ j Rf{uj{ip),0)\sinip\dip = ^ j Rf{uj, 0 )\uj-e 2 \duj, 

0 0 si 


which is the equation ffTT]) for n = 2. 

In order to prove the proposition for n > 3, we need two auxiliary results. 

Lemma 23. For ipo G (0, vr/2), ip G (0, tt), and n > 3, we define 


^(V’o,V’) 


(cos^ Ipo — COS^ -0)^” 
{sinip)^~^ 


Then, for any f G S{MT), 


(31) 


TT-'ljjO 

f (COSf 

/ Cf{0,en,'ip)g{'ipo,'4’)dip = —— / Rf{{cosipo)u,sinipo),0)du. 

IpO 

(32) 

Proof. The idea of the proof is to exhaust the exterior volume of two opposite 
cones having a common vertex in two ways. The hrst is by taking a family 
of vertical cones whose vertices are at the origin and opening angles vary 
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from -00 to vr — '^ 0 - The second is to consider a family of hyperplanes passing 
through origin and are tangent to the vertical cone having vertex at the origin 
and opening angle tfjQ (See Fig. [6]). 


Z 



Let the functions / and g be given as in the lemma. We can split the 
integral on the left hand side of equation fl5^ into two parts to get 


TT —-00 

[ f C/(0, e„, ^)# 


ho 


ho 

TT—ho 


+ j Cf{0,en,'ip)g{'il^o,'ip)d'4!. (33) 

7r/2 


By the dehnition of the vertical cone transform (jl]), for the hrst term on 
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the right hand side, then 


t/2 




^0 


7r/2 oo 



/(p(sin'0)a;, pcos'0)(psin'0)”' 'ifj)dudpd'ip. 


i/'o 0 5"-2 


If we make a change of variables in the integral with respect to p by 
letting z = pcos'ip, we have 


tt/2 


C/(O,en,' 0 )fi'(' 0 o,' 0 )# 


ho 


tt/2 



dz 


/(ztamjjcu, z)(ztamjj)'^ g(ijjo, tjj^duj - -d^jj. 

cos'll! 

ho 0 

Now, if we let r = ziamp, then dr = z sec^ ipdip^ and since 

sec^ -0 — sec^ 0o tan^ 0 — tan^ 0o ~ f 0o 


cos^ 00 — cos^ 0 = 


we have 5 f( 0 o, 0(r, z)) = 
Thus, 


sec^ 00 sec^ 0 sec^ 0o sec^ 0 sec^ 0o sec^ 0 ’ 

tan^ 0o)*'’^“"‘^/^ 


(r sec 00 ) 


n—4 


7r/2 


C/(O,e„,0)5((0o,0)d0 


ho 


OO OO 


= (cos0o) 


n—4 



J fz{ru){r'^ — z'^ tan^ ^l!oY"^ 

0 5^—2 2 : tan 'i/’o 

Then, using the identity fl2^ . we obtain the following relation between 
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the cone transform of / and {n — l)-dimensional Radon transform of fz- 

7r/2 

t rns /" f 

Rfzioj, —ztamjjo)dudz 


fcos 


do 


1 ^" 



(cos V’o) 

1 


n—4 


n—31 



0 5^-2 

fz{x)5{x ■ u + ztaia'ipQ)dxdudz. 


0 

- 1 ; 


Now, since (5 (A(m — a)) = A — a), we get 


t/2 


Cf{0,en,^p)g{^|Jo,^p)d^p 


do 


(cos' 0 o)’^ ^ 
1 5^-3 1 



f{x,z)6{x- (cos'ipo)u + z sm'ipQ)dxdzdu. (34) 


5"-2 0 R"-l 


For the second term of the right hand side of (13^ , we change the variable 
V’ by TT — -0 to get 

TT—ipo 7r/2 

/ Cf{0,en,'ip)g{'ipo,'4’)d'ilJ= / C/(0, e,,, vr - ?/>)5((V'o, tJ" --0)# 


7r/2 


do 


7r/2 



/(p(sin'0)a;, —pcos'0)(psin'0)"' ‘^g{'ijJo, il))dojdpd'il). 

V’o 0 s^~‘^ 

Again we change variables hrst by letting z = p cos ip and then r = z tan ip 
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to obtain 


■n—ipQ 


C'/(0,e„, 


t/2 


oo oo 


= (cos?/>o) 


n—4 



f-z (ru) (r^ — z'^ tan^ '^o) ^^^^rdrdudz 


0 5 ^“^-2 tan-i/tQ 


(cos'0o)"' 


n—31 



Rf_z{^, zta3i'ipQ)dudz, 


0 


where the last equality follows from the identity (12^ . Again, by the definition 
of the Radon transform, and 5{X{u — a)) = \~^6{u — a), we get 


-R/_^(a;, z tan'ipQ)dudz 


0 S’i-2 



= cos-00 



f{x, —z)S(x ■ (cosi/^o)^ ~ zsinil^ojdxdzdcj 


gn-2 0 R"-l 
0 


= cos-00 


f{x, z)S{x ■ (cos + 2; sin %l)Q)dxdzdu. 


Sn-2 _ooRn-l 


Thus, 


TT—lpO 


Cf{0,en,'ip)g{'ilJo,'ijj)d'4! 


7r/2 


(cos'^o)"' ^ 
I ^^-31 


f{x, z)6{x ■ {cos'ipo)u + zsm'ijjQ)dxdzdu. 


gn-2 _ooR^-l 


(35) 


Now, using fl34l) and fl3^ for the first and second terms in the equation 
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dMl), we obtain 


TT — lpO 


Cf{0,en,'4’)g{^po,'ilj)d'ijj 


V’o 


(cosV’o) 

1 ^' 


n—3 


n—3 1 



f{x)6{x- {{cos'ipo)^,s^^4’o))dxdu. 


Finally, observing that 


f{x)6{x ■ {{cos'ipo)(^,sm'ipo))dx = Rf {{{cos'ipo)(^,sm'ipo),0), 


we have 
TT—yo 


Cf{0,en,^/j)g{i/Jo,i^)d^/j = — [ Rf {{cosilJo)uj, sin7po),0)du. 


15 ^- 31 

tpo 5’^“2 

Hence, we get the result. □ 

Lemma 24. Assume that n> 3. Let g{'ipo,ip) be given as in (l3T|l and define 

(cos^'i/'o — 


h{'i/jo,fi) = 


(sin fi) 


n—3 


Then, 


TT — lpO 


dfio 


ho 


Cf{0,en,'fi)h{fio,fi)dfi 

TT-ho 

= (2 - n) COS-00 sin-00 j Cf{0,en,fi)g{^|Jo,^|J)d^jJ. 

ho 


dh 


Proof. As 7 -— {fio, fi) = (2—n) cos V’o Amjjog{fiQ, fi), utilizing Leibniz integral 
ofio 

rule and noticing that h{'fio, tt — ipo) = h{'ipQ, iPq) = 0 gives the result. □ 
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Proof of Proposition\^ n > 3. By Lemmas 12^ and 12^ we have 

TT — lpO 

Cf{0,en,fj)h{fjo,7p)dfj 


d 


d'lpo 


■00 


TT-tpo 


= (2 - n) COS V'o sin ^/>o / Cf{0, en,fj)gifJo,fj)dxp 

-R/(((cos' 0 o)i^, sin fjo), 0)du. 


00 


(2 — 'n)sm'0o(cos'0o)"' ^ 


1 


n—31 


Sn-2 


Integrating both sides with respect to fjo from 0 to 7 r/ 2 , we obtain 

TT 

[Cf {0, Bn, fj) sin'ljJd^jJ 


n — 2 
I ^>^-3 1 

n — 2 
I ^>^-3 1 


7r/2 



i?/(((cos^/>o)i^,sinV’o),0)(ia;sin^/>o(cos^/>o)" #o (36) 


0 

TV f2 



-R/(((sin 0 )a;, cos 0 ), 0 ) cos 0 (sin 0 )"' ^dudf), 


0 


71 


where we changed the variable by letting 0 =- t/jq. On the other hand, 

TT 

letting 0 = 00 + —, we have 


0 /( 0 , e„, 0 ) sin 0 (i 0 

7r/2 


n — 2 


n—3 1 


1 5^ 


n — 2 



Rf {{{cos 7jjo)u, sin 00 ), O)(ia;sin0o(cos0o)"^ ^d0o 


0 5»j-2 


1 5^ 


'n—31 



-R/(((sin 0 )a;, — cos 0 ), 0 )(— cos 0 )(sin 0 )"^ ^dudcp. 


1x12 5" 


27 








Now, due to evenness of Radon transform, we have 

Rf {{{sin(j))uj, — COS0), 0) = Rf{{{— sin0)(—cn), — cos0), 0) 

= R/(-((sin0)(-a;),cos0),O) = Rf {{{sin (j)){-uj), cos 

Since the Lebesgue measure is rotation invariant, we obtain 


TT 

j Cf{0,en,^p)sin^jJd^p 

0 


n — 2 


\S^' 


n—31 



Rf {{{sin (j))u, COS0), 0)(— cos0)(sin(/))” ^doodf. 


7r/2 S"-2 


(37) 


Summing (l36|) and (1^ . we conclude that 


j Cf{0,en,'f)simpd'ijj 
0 

TT 

= 2|gn-^3| J J ^/(((sm0)a;,cos0),O)|cos0|(sin0)'‘"2^a;(i0 

0 

n — 2 C 

= 2|gn-3| J Rfi(^,0)\a-en\da. 

5n-l 

Finally, application of the formula flTTll and r( 2 ; + 1) = zV{z) gives the 
result. □ 


8 Conclusions and Remarks 

In this paper, various relations between the general (overdetermined) cone 
transform and Radon and cosine transforms and spherical harmonic expan¬ 
sions are explored. Several inversion formulas for the cone transform are 
obtained, some of which of hltered backprojection nature. Examples of re¬ 
constructions from synthetic Compton camera data are provided. 

Some additional remarks: 
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• In order not to distract from the main point, the source intensity dis¬ 
tribution function / is assumed to be of the Schwartz class, 

In fact, the cone transform of / is well-dehned even when we assume 
integrability of / on each cone. The formulas obtained here can be 
extended by continuity to much larger function spaces. For instance, 
for the inversion formula flT^ to hold, it is sufficient that the function 
Cf is {n — l)-times differentiable with respect to w, and to this end it 
suffices to assume the function / be (n — l)-times differentiable. As a 
condition of decaying, assuming that f{x) = 0{\x\~^) for some N > n, 
is sufficient. 

• Although we do not explicitly present the adjoint of the cone trans¬ 
form, both Theorem |2] and Theorem [8] provide hltered back projection 
type inversion formulas for the cone transform as, in both cases, we 
recover the function at a point u using a weighted averaging of its cone 
transform over cones having vertex at u. 

• Let us address the comparison of inversion formulas of Theorems [2] and 
IHl Both of them involve integrating the data with respect to 'ijj and 
[3 and hltration by the same Riesz potential. The difference is that 
in Theorem |2] the measure of integration is fi{(3)d‘ipdf3 with arbitrary 
function /i of mass 1 (e.g., a (5-function), while the formula of Theorem 
[8] holds only for the measure sm{‘ip)d‘ipd(3. 
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